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Abstract 


There is strong interaction between heat/mass transfer and electrode dynamics in an air-breathing proton exchange membrane fuel cell (PEMFC). 
To investigate the heat/mass transfer characteristics of PEMFC, their effect on oxygen transport and then on the performance of the cell, a 
coupled three-dimensional (3D) mathematical model has been developed for an air-breathing PEMFC using non-dimensional heat/mass transfer 
coefficients, and a detailed interpretation of heat/mass transfer parameters in it has been done. The full elliptic Navier Stokes and energy equations 
are simultaneously solved in the composite domain with commercial CFD tool STAR-CD based on the finite volume numerical method. And the 
model has been confirmed by experimental results. From the numerical simulation results obtained, a simple equation has been suggested on the 
function of the sum of Grashof number and diffusional Grashof number for the calculation of dimensionless mass transfer coefficient Sherwood 
number. The equation is believed to be useful for a profound understanding of the nature of oxygen transport limitation coupled with the heat/mass 
transfer in the air-breathing PEMFC, and further for the optimum design of the air-breathing PEMFCs. The concentration over-potential as a 
function of limited oxygen mass transfer rate, as well as profiles of velocity, temperature, reactant and water concentration, current density are 


presented and discussed to analyze the coupling problems between heat/mass transfer and electrode dynamics in an air-breathing PEMFC. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Due to high electrical efficiency, flexibility with respect to 
power and capacity, long lifetime and good ecological balance, 
fuel cells have become one of the most interesting alternatives 
for clean power production in automotive applications and for 
distributed power generation. In addition to fuel cells for station- 
ary or automotive applications, fuel cells also have the potential 
to complement or to substitute batteries in portable applications. 
The air-breathing PEMFCs operating at an ambient temperature 
without humidifier and air compressor are one of the most likely 
candidates for miniature fuel cells. A typical structure and prin- 
ciple of an air-breathing PEMFC that is used in our laboratory 
are illustrated in Fig. 1. 

In this single air-breathing PEMFC, the structure consists of 
seven components as shown in Fig. 1. The electrolyte used here 
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is Nafion 112, and the catalyst is platinum supported carbon. 
With two electrodes, the electrolyte and catalyst are assembled 
into a sandwich structure to form membrane-electrode-assembly 
(MEA). The cathode bipolar plate is machined to straight chan- 
nel and anode plate designed as serpentine groove as shown in 
Fig. 1. 

It has been known that the oxygen transport limitation plays a 
great role in the performance of air-breathing fuel cells from lit- 
eratures [1-3]. In the air-breathing PEMFC, air flow into channel 
is driven by the effect of buoyancy, which is distributed through- 
out the air and is associated with the tendency of air to expand 
when heated. Air is heated due to irreversible losses when the 
electrochemical reaction occurs in the catalyst interface. The 
heated air becomes less dense and flows upwards. This process 
is called free convection. Thus the temperature gradient and 
electrochemical reaction has a remarkable impact on air flow 
rate. In addition, the liquid water also has a significant impact 
on the oxygen transport by blocking the gas way in the chan- 
nel and backing layer when the air-breathing PEMFC works 
in high ambient humidity. The cell performance is decreased 
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Nomenclature 

a constant of the binary diffusion coefficient 

Acv specific area of control volume (m7!) 

b exponent of the binary diffusion coefficient 

Cò, oxygen concentration in ambient air (kg m~?) 

D mass diffusion coefficient (cm? s7!) 

F Faraday’s constant (96,487 C mol`!) 

g gravity (9.8 m s7?) 

I local current density (A cm~?) 

Tim concentration limiting current density (A cm~?) 

P exchange current density (A cm7?) 

K thermal conductivity (J s7! m7! °C!) 

Kp the permeability of the porous medium (m~?) 

l characteristic length (m) 

M molecular weight of gases (kg mol™!) 

mo, oxygen mass flux in cathode catalyst surface 
(kg m~? s7!) 

na electro-osmotic drag coefficient 

P pressure (Pa) 

Po critical pressure (atm) 

R gas constant (8.314 J mol -1 K7!) 

S source term in governing equation 

T temperature (K) 

Ter critical temperature (K) 

tp diffusion layer thickness (mm) 

fem membrane thickness (mm) 

u fluid velocity vector (m s7!) 

Ve cell voltage (V) 

Vintrinsic intrinsic average velocity (m s7!) 

Voc open circuit voltage (V) 

w mass fraction 


Greek letters 


net water molecules per proton flux 
coefficient of thermal expansion i) 
porosity for a porous medium 
over-potential (V) 

dynamic viscosity (kg m7! s7!) 
change in... 

gas phase density (kg m~*) 


Superscripts, subscripts 


anode 
activation 
cathode 
hydrogen 
water 

gas mixture 
membrane 
ohmic losses 
oxygen 

wall 

initial condition 


7 


Mi = 


b. Air channels c. Fuel channels 


Outlet air reservoir 


(1) Bipolar plate 
(2) anode GDE 
(3) anode catalyst layer 


(4) membrane 
(5) cathode GDE 
(6) cathode catalyst layer 


Inlet air reservoir 


PE d. Numerical simulation domain 
a. Structure and principle 


Fig. 1. Scheme of geometry and calculation domain of an air-breathing PEMFC. 


extraordinarily when the oxygen partial pressure is insufficient 
in catalyst layer. Therefore, the performance of an air-breathing 
PEMFC is greatly influenced by heat/mass transfer as well as air 
free convection and water transport. And all these processes are 
very closely coupled. Understanding these transport processes 
and the effects of them on the performance of electrochemi- 
cal reaction in an air-breathing PEMFC is evidently important 
for optimization of the fuel cell design and operating condi- 
tions. Hence a comprehensive modeling of these processes is an 
effective way of doing optimization. 

Many models focused on general forced convection PEMFC, 
most of them are one-dimensional (1D) flows, Springer and his 
co-workers developed an isothermal, one-dimensional, diffusion 
model. At the same time, Bernardi and Verbruge [4] developed 
a heat and water management model in dimensions normal to 
the reactive interface. Nguyen and Springer [5,6] developed 
a one-dimensional mathematical model in dimensions along 
the channel. One-dimensional models restrain the application 
to simple geometry and isothermal assumption even though 
they have provided useful insight and fundamental understand- 
ing. To predict the performance of fuel cells, to investigate the 
operation of fuel cells and to optimize the distribution of gas, 
two-dimensional and three-dimension models were developed 
recently. Hsing and Futerko [7] developed a two-dimensional 
model of coupled fluid flow, mass transport and electrochem- 
istry of PEMFC operating without external humidification of 
hydrogen and air streams. Dutta [8,9] described a steady state, 
isothermal, three-dimensional and single phase PEMFC model 
used to predict the performance of a fuel cell. In addition, 
recently there are several works about heat/mass transfer, and 
air-breathing PEMFCs have been published. Jeong et al. [10] 
investigated the characteristics and optimized the cathode struc- 
ture of air-breathing PEMFCs through experiment. And the 
influence of gas diffusion layer properties and cathode catalyst 
loading on the cell performance are examined in their works. 
Further, in their another literature [11], they investigated exper- 
imentally the effects of cathode open area and relative humidity 
on the air-breathing fuel cell performance, and found that the 
single cell with the open area of 77% exhibited the highest 
performance. Ramousse et al. [12] presented a fuel cell model 
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in the condition of forced convection considering gas diffusion 
in porous electrodes, water diffusion, electro-osmotic transport 
through the polymeric membrane and heat transfer in both the 
membrane electrodes assembly (MEA) and bipolar plates. The 
results showed that thermal gradients in the MEA could lead 
to thermal stresses at high current densities. Arato and Costa 
[13] discussed the interaction between diffusive and forced 
flows expressed in terms of Peclet numbers and the overall dif- 
fusive resistance in terms of Sherwood numbers in a forced 
convection PEMFCs with emphasis on the oxygen diffusion 
limitation. The results of this approach were then used in the 
discussion of the performance of different geometric arrange- 
ments, such as traditional, interdigitated and serpentine cells. 
Shan and Choe [14] presented a dynamic PEMFCs stack model 
taking into account heat transfer, and analyzed start-up behaviors 
and the performance of the stack. Wang et al. [15] developed a 
two-dimensional numerical model for an air-breathing PEMFC 
including water transportation and electrochemical reaction in 
catalyst layer. In this model they discussed the effect of cathode 
open size on cell performance and suggested a best opening for 
the cathode channel. 

In spite of the fact that the above-cited references have 
stressed the modeling and calculation for fuel cells to analyze 
the various phenomena, none of them addresses the full three- 
dimensional solution for heat/mass transfer, fluid dynamics and 
electrochemistry in an air-breathing PEMFC. It is extremely 
challenging to derive and to solve those kinds of strongly cou- 
pled model equations. In previous studies, Wang et al. [16,17] 
have investigated the effect of cathode plate configuration on cell 
performance by developing a three dimensional CFD model for 
an air-breathing PEMFC. Experimental data has also been col- 
lected to confirm the model result and the cell performance under 
different air humidity, the temperature on the carbon surface. 

In this work, we developed a model that could analyze in 
detail the important characteristic of heat and mass transfer 
in an air-breathing PEMFC. A correlation between Sherwood 
number and Grashof number has been derived on the basis of 
the numerical results, and the correlation could greatly sim- 
plify the calculation of the limiting current density without the 
development and solution of a full N-S equation. 

The computation domain comprised an actual physical 
domain between two bipolar plate, as well as two relatively large 
domains placed on the upstream of channel entrance and on the 
downstream of the channel exit as shown in Fig. 1. These I- 
type composite domains resemble open, unbounded spaces that 
serve to accommodate the diffusion phenomena by momentum 
and energy that could occur outside of the channel. Using com- 
putational fluid dynamics technique, the conservation equations 
have been solved using velocity, temperature, concentration of 
gas component, over-potential and current density fields to ana- 
lyze the effect of heat and mass transfer on the performance of an 
air-breathing PEMFC with the perspective of its optimization. 


2. Model descriptions 


The air-breathing PEM fuel cell model presented here is a 
comprehensive three-dimensional, non-isothermal, steady state 


model for explaining the complicated process of heat and mass 
transfer, electrochemistry and fluid dynamics in an air-breathing 
PEMFC. 


2.1. Model assumptions 


A complete air-breathing PEMFC is an extremely complex 
system. Therefore, to obtain the solution of three-dimensional 
model of an air-breathing PEMFC it is necessary to do some 
reasonable simplifying assumptions. In this model following 
assumptions are used: 


e The water exists in the fuel cell as vapor water. It is reasonable 
because the water partial pressure is not high in air-breathing 
PEM fuel cell without external humidification but with low 
ambient air humidity. 

e Ohmic heating and ohmic resistance in the bipolar plate and 
in the carbon gas diffusion layer are neglected due to high 
thermal conductivity and electrical conductivity. 

e Membrane is assumed to be solid for faster numerical con- 
vergence. 


2.2. Governing equations 


In the gas flow field, the flow is a laminar steady state flow. 
The governing equations of the fluid flow represent mathemati- 
cal statements of the conservation laws of physics. Considering 
the effect of porous medium, from the Dupuit—Forchherimer 
relationship (u = €Vintrinsic) [18], the fluid conservation equa- 
tions in porous medium can be derived as following. The detailed 
arguments can also be found in Wang’s work [19]. 


2.2.1. Continuity equation 
The mass conservation can be stated for steady state flow as 
following: 


V. e(pu) = Smass (1) 


€ is porosity for a porous medium, which is zero in channel 
and Sass is the mass source or sink because of electrochemical 
reaction, which is sum of source terms for oxygen, nitrogen, 
vapor water in cathode and for hydrogen, vapor water in anode. 


2.2.2. Momentum conservation 

Following momentum equation is derived from Newton’s 
second law: 

ae z euu 

V . (puu) = —VP + V - (uVu) + pg — (5) (2) 
Here the terms from left to right are, respectively, convective 
terms, pressure drop terms, diffusive terms, buoyancy effect 
terms according to Boussinesq approximation and terms for 
pressure drop due to small permeability of porous medium based 
on Dancy’s law. Kp is the permeability of the porous medium. 
The dynamic viscosity of gas mixture can be calculated from 
following formula [20]: 
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(3) 


where Y; and Y; are molar fraction of gas components i and j, pi 
is dynamic viscosity of gas component i, i and j represent either 
oxygen, nitrogen or vapor in cathode, and hydrogen or vapor 
in anode, the non-dimensional parameter g;; can be calculated 
from following formula [20]: 


1 M,\ 712 wi \ 2 / M" f 
J J i 


where M; and Mj; are molecular weight of gas components i and 
j. 


2.2.3. Energy conservation 
The energy equation is derived from the first law of thermo- 
dynamics: 


V- e (pcp) = V : (KVT) + Sr (5) 


Here K°f is the effective thermal conductivity of gas mixture 
in porous medium, which can be yielded as a simple function 
of the gas mixture thermal conductivity in nonporous system. 
This approximate equation is derived from the analogue with 
the calculation of the effective diffusion coefficient to simplify 
the numerical process. 


Kg = Kel (6) 


cp is the specific heat capacity at constant pressure and K is 
the thermal conductivity of gas mixture, which can be estimated 
with a function of the individual gas component’s thermal con— 


10.13 x a(T/ s/Ter,iTer,j) (Peri Per, j)/3 (Ter,iTer, Q/M) + (1/M p)? 


The momentum equation and energy equation are coupled 
via equation of state of ideal mixtures. 
PMm 


p= RT (10) 


where Mm is the molecular weight of the gas mixture which 
can be calculated from mass fraction of gas component in gas 
mixture. 


1 
~ i(wi/ Mi) 


where w; is the mass fraction of gas component and M; is the 
molecular weight of gas component i. 


Mm (1) 


2.2.4. Species mass transport 
The mass fraction of gas component can be calculated from 
following transport equation: 


V - (eowjit) = V - (D$ËV pw) + S; (12) 


Here pe is effective molecular diffusion coefficient of gas com- 
ponent i in the porous medium and have a relationship with 
gas diffusion coefficient in nonporous material D; as stated in 
literature [18] 


pe = Deg? (13) 
D; can be calculated from binary diffusion coefficients. 
_ A/M = wi) 

> j¢i(wj/DijMj) 
where the binary diffusion coefficient Dj for gas components i 
and j is calculated from Slattery and Bird estimation [20] 


(14) 


i 


ij 


ductivity weighted by the mass fraction from following formula 
according to the methodology of STAR-CD: 


N 
=X wigi (7) 
i=l 


¢ represents the specific heat capacity or the thermal conductiv- 
ity of gas mixture, ¢; represents the specific heat capacity or the 
thermal conductivity of gas component i. 

Sr is energy source term which represents the rate of increase 
or decrease of energy due to heat generations or consumptions. 
In the catalyst layer, because of the heat transformed from irre- 
versible over-potential, the source terms in cathode and anode 
are, respectively: 


St = Acvina, Sp = Acv!nc (8) 
In membrane, heat generated from ohmic resistance. 
srem = Acv ohm (9) 


Here Acy is specific surface area of the control volume. 


15 
Pp (5) 
The value of parameters used in above equation can be found 
in Table 1. According to Faraday’s laws, the source term S; in 
species mass transport equation for oxygen and hydrogen is: 
—IM, Acv 
S.= 16 
= (16) 
where r is oxygen or hydrogen, M, the molecular mass weight 
of reactant r, n is 2 for hydrogen and 4 for oxygen. 
Vapor water source term can be calculated from the balance 
of water transport. Vapor water source terms in anode and in 
cathode are: 


_ MmoAcva 


a = 
Sho = - ee" (17) 
Table 1 
The value of parameters used in the equation for binary diffusion coefficient 
Species A Ma AAB bap Poa Tea Species B 
O2 31.999 3.64E—4 2.334 49.7 154.4 H20 
N2 28.013 2.75E—4 1.823 33.5 126.2 Or 
H20 18.015 3.64E—4 2.334 221.2 6473 M 
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se _ IMmoAcv aMy,0Acv 
H20 2F F 


(18) 


Here net water molecules per proton flux œ can be yielded from 
following water transport equation through membrane: 


a= na FDy(Cw,c = Cw,a) (19) 


Itmem 


The first term in the right of the equation is water molecules 
from anode to cathode by electro-osmotic drag caused by proton 
transport. The second is the diffusion term of water molecules by 
concentration difference between anode and cathode side, based 
on a single-step linear water concentration difference assump- 
tion. Cw, and Cw, are, respectively, water concentrations in 
cathode and anode. The electro-osmotic drag coefficient ng, and 
the water diffusion coefficient in membrane Dw can be obtained 
from the function of the water content in membrane. The detailed 
membrane model has been presented in literature [15]. 
The mass fraction of nitrogen can be obtained by 


wn, = 1 — wo, — wHo (20) 


In above equation, / is the current density that can be calculated 
from following equations [6,8]: 


Ve = Voc — Nact — Nohm — conc (21) 

< l (22) 
a = —— In 

Nact 05F 70 Po, 
t I 

Nohm = a (23) 
Omem 

—4 I 
Neonc = —3.4 x 10°-"T In| 1 — — (24) 
lim 


Here the membrane conductivity om can be obtained as a func- 
tion of temperature and water content in membrane [9], the 
concentration over-potential nconc is related to the oxygen trans- 
fer rate in the cathode channel. And then the limiting current 
density [jim can be calculated from oxygen transport limitation 
according to literature [21]. 


4FDo Ce. Sh 
Kim = ———— (25) 
tD 
where the dimensionless mass transfer coefficient Sherwood 
number (Sh) is calculated by following equation based on the 


definition in literature [20]: 


= (tp/Do,)mo, 


Sh = 
p(we, = Wo.) 


(26) 


2.3. Numerical scheme 


The calculation procedure is carried out in commercial CFD 
tool STAR-CD with the parameters, boundary conditions and 
geometry dimension listed in Table 2. The initial source terms 
for governing equations are obtained from average current den- 
sity and calculated over-potential at first sweep of iterations. 
The governing equations are solved based on the SIMPLE algo- 
rithm and the monotone advection and reconstruction scheme 
(MARS) using finite-volume approach. The partial pressure of 
gas component and the oxygen mass transfer coefficient are 
obtained from the solution and are used to calculate the current 
density, over-potential and the net water molecules per proton 
flux by user-defined subroutine based on the finite-difference 
discretization technique. Then with the new source terms calcu- 
lated using newly obtained values the governing equations are 
solved iteratively until the convergence. 


3. Results and discussions 


In the following paragraph the results are presented for an air- 
breathing PEMFC operated in the condition shown in Table 2, 
with an electrode area of 34.4 cm?, the anode and cathode chan- 
nel arranged as intercross flow mode. The analysis has been 
made to investigate the effects of heat and mass transfer, water 
profile on the performance of the air-breathing PEMFC. 


3.1. Model validation 


The fuel cell terminal voltage was calculated for series 
of cell current density based on the operating conditions 
shown in Table 2. The experiment has been carried out in our 
laboratory with the same cell geometry specified to validate 
the calculated voltage values. Fig. 2 compares the calculated 
voltage for this air-breathing PEMFC as a function of current 
density with the experiment results. Generally the computed 
voltages for several current densities shown in Fig. 2 are in good 
agreement with the experimental data. However, the calculated 
performance is little higher than the experimental result. This 
discrepancy is believed to be due to the experimental errors 


Table 2 

Values for parameters and boundary condition used in the model 

Parameters Symbol Value Parameters Symbol Value 

Porosity in carbon paper and E 0.4, 0.28 Cathode channel height and width he, We 0.3 cm 
catalyst layer 

Average current density Taye 0.15 A cm~? Relative humidity in air E 63% RH 

Ambient temperature, pressure Tp, P 300 K, 101,325 Pa Membrane thickness tm 0.0508 mm 

Anode channel width, height Wa, Na 0.1 cm Thickness of Diffusion layer, catalyst layer tats 0.026 mm, 0.001 mm 

Hh inlet temperature Ta 299K Oxygen/nitrogen ratio in air ¢ 0.21/0.79 

Hp flow rate in inlet Np, 7.198E—08 kg s7! Cathode channel length Le 47mm 

Permeability in GDL Kp 1.8E—18 m? Permeability in catalyst layer Kp 1.76E—12 m? 
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Fig. 2. Comparison of polarization curves between the numerical prediction and 
experimental results. 


caused by difficulties in controlling of ambient humidity and 
temperature. 


3.2. Effect of heat transfer 


To understand the effect of the heat transfer, the calculated 
profiles of temperature, Grashof number (Grr), velocity and 
current density are presented here at the current density of 
0.15Acm7!, which is estimated to be an optimal operation 
current load from Fig. 2. 

Fig. 3 displays the temperature profile in the MEA region. 
Because of low thermal conductivity in the diffusion layer and a 
low heat transfer coefficient in catalyst layer, the gas temperature 
in the catalyst layer is the highest. Given the temperatures in the 
cathode and anode diffusion layer in the figure, those tempera- 
ture profiles in this air-breathing PEMFC is different from those 
of general forced convection PEMFC where the temperature of 
the anode is a little higher than in cathode. This is apparently 
due to the fact that in an air breathing PEMFC, the cathode, 
which has an air channel, can release more heat to the ambi- 
ent air than the anode. Observing the temperature distribution 
in the catalyst interface, the temperature is high in the center 
due to the presence of thermal diffusion and convection in the 
air channel. And the temperature is high in the region near the 
anode inlet and outlet because there is a strong electrochemical 
reaction caused by affluent hydrogen or water. The region near 


39.76 
39.61 
39.45 
39.30 
39.15 
34.99 
34.64 
34.68 
34.53 
34.38 
34.22 
34.07 
33.92 
33.76 
33.61 


Fig. 3. Temperature distribution in different layers of MEA (°C). 


the outlet of the anode should be rich of water because of con- 
tinuous water generation from electrochemical reaction. For this 
air-breathing PEMFC operated under low humidity condition, 
higher water concentration could be better for electrochemical 
reaction. Fig. 4 shows the temperature profile of the cathode 
channel. The thermal boundary layer can be clearly observed. 
In the air-breathing PEMFC, air is heated from the bottom near 
the MEA surface, the heat transfer includes the heat transport 
from the MEA to the air channel and to the solid plate by con- 
vection and conduction, and the heat transport between air in the 
channel and solid plate. To simplify the model the heat transfer 
resistance between different regions was not taken into account. 
Evidently there should be a trade-off between the heat trans- 
fer to solid plate and to air in the channel in order to improve 
the free convection in the channel. But it is difficult to give a 
clear definition to the optimal ratio for the heat transfer rate 
between to air and to solid plate due to coupled effects between 
these two. However, in the air-breathing PEMFC, what is more 
important than the heat transfer itself is the mass transfer rate 
that is limiting the performance of the cell. The heat transfer put 
an important effect on the mass transfer through the buoyancy 
forces. The dimensionless parameter Grashof number defined 


Fig. 4. Temperature distribution in cathode channel (°C) (=-= is the view direction for the section plot). 
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Fig. 5. The distribution of local Grashof number in cathode catalyst surface. 


as the strength of buoyancy effect due to temperature differ- 
ence could greatly contribute to the mass transfer rate in the air 
breathing PEMFC. Fig. 5 shows the distribution of local Grashof 
number. According to Bird et al. [20], the Grashof number can 
be calculated from a simple formula, 


_ ghrp? (Tw — Tol? 
u2 

The distribution of the Grr is not uniform due to uneven elec- 
trochemical reaction. Comparing Figs. 3-5, it can be observed 
that the heat transfer takes place almost throughout the full length 
of the channel. And in the middle of plate, in the opening of the 
furthest right and left channel near the anode inlet and outlet, 
the heat transfer has a significant effect on free convection of the 
air. This also can be proved from the profile of velocity shown in 
Fig. 6. The center channels have a longer entrance length than 
others. The flow is not developed yet and shows strong mass 
and heat transfer in the entrance region. In the middle channel 
marked by ‘3’ in the figure has longest entrance length. In the 
channels, marked ‘1’, ‘2’ in the figure, air flows into channel 
not only from the bottom inlet but also from the top opening 
due to characteristic of heat transfer and electrochemical reac- 
tion. Fig. 7 shows the current density distribution. Similarly the 
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Fig. 6. The velocity profile in cathode channel, ms~!. 
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Fig. 7. The current density distribution in catalyst surface, Acm~?. 


Grr (27) 


electrochemical reaction is increased in the channels ‘1’ and ‘2’, 
and the regions near the opening of the cathode channel. From 
Figs. 3-7, it can be observed that there exsits strong interaction 
between heat transfer and electrochemical reaction. Effective 
heat transfer in the air channel can improve the performance of 
air-breathing PEMFC by causing uniform gas distribution, and 
by sustaining a good heat transfer rate between solid and fluid, 
and observed in the long thermal and velocity entrance length 
in the cathode channel. 


3.3. Effect of mass transfer 


The performance of air-breathing PEMFCs is strongly 
affected by free convection in air channel. When air supply to 
channel is not sufficient, the oxygen partial pressure in the cata- 
lyst layer is accordingly low; the performance of cell decreased 
by the function of Nernst equation, this phenomenon is called 
mass transport limitation losses or concentration losses. If the 
oxygen concentration is less than some level in catalyst layer, 
the cell will stop working completely. Fig. 8 is oxygen mass 
fraction distribution in cathode side. It shows that the oxygen 
mass fraction is much lower than in ambient air even though 
the cell works in such low current load. The calculated air mass 
flux to catalyst layer by convection is 6.818E—7kgs~!, and 
hence approximately, 1.32E—7kgs~! oxygen flows into cata- 
lyst layer by convection considering the oxygen mass fraction 
in air. Moreover, there is a little oxygen diffusion to catalyst 
layer. Total oxygen supply should sustain the electrochemical 
reaction at this average current density. The air-breathing PEM- 
FCs can not work at high current density because of limited 
oxygen supplement. The current density at which the oxygen is 
used up is called limiting current density where the cell cannot 
work beyond. In the air-breathing PEMFC, the limiting current 


Fig. 8. Oxygen mass fraction in cathode side (G.D.L. is abbreviation of diffusion 
layer, R.L. is abbreviation of catalyst layer). 
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density is closely associated with mass transfer rate to catalyst 
layer. The mass transfer in the air-breathing PEMFC is decided 
not only by free convection related to heat transfer, but also by 
the oxygen diffusion caused by concentration gradients at the 
electrochemical reaction zone. Similar to the heat transfer, the 
diffusional Grashof number due to concentration difference can 
be derived as following for mass transfer [20]: 


2 3 
w —w l 
Cig = gyp~( On. 02,0) (28) 
u 


As thermal coefficient of volume expansion in free convection 
heat transfer, mass expansion coefficient is 


__(1)(_¥ 
r=- (5) (sco) o 


For mass transfer between three ideal gases, the expression for 
them can be further derived as 


= My,0Mn, — My,0Mo, 
My,0Mn, Wo,,w + Mo, MN, WHO, w 
+ Mo, My,0(1 — wH,o,w. — Wo2,w) 


Yoz (30) 


Fig. 9 shows the distribution of diffusional Grashof number. 
On the contrary to heat transfer, in the right and left furthest 
channel, the diffusional Grashof numbers are low due to small 
concentration differences. Therefore, considering the heat trans- 
fer and mass transfer simultaneously, it can be found that the 
mass transfer rate is affected by the combination of heat trans- 
fer and mass transfer. Considering the heat transfer and mass 
transfer by the free convection in the channel and in the porous 
medium, there are few equations applicable till now especially 
for mass transfer rate in this case. Literature [20] derives it from 
analogy of heat transfer coefficient. However, more reasonable 
values could be derived from the air flow rate which can be 
calculated from solution of full Navier-Stokes equation. 

Fig. 10 shows the Sherwood number distribution, which 
represents the oxygen transport characteristic in cathode. Com- 
paring with current density distribution in Fig. 7, it is found 
that these two figures have very similar pattern. Therefore, it 
can be concluded that the mass transport plays major role in 
cell performance. Comparing the profiles for Grashof number 
and diffusional Grashof number with those of Sherwood number 
distribution, it can evidently be found that the free convection 
due to temperature difference and the diffusion due to concen- 
tration gradients has great impact on the oxygen transport, and 
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Fig. 9. The distribution of diffusional Grashof number in cathode catalyst sur- 
face. 
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Fig. 10. Sherwood number distribution in cathode catalyst surface. 


there are some coupled relationships between the Sh and Grm, 
Grr. Therefore, it is more desirable that Sh is calculated by 
the solution of full N-S equation, however, the complexity and 
difficulties encountered to solve the equation limited the appli- 
cation of this method. Therefore, it is necessary to develop a 
simple equation to calculate the Sh. A correlation of Sh and 
Grm, Grr has been suggested for the vapor natural convection 
model including heat and mass transfer simultaneously [20]. 


Sh = Ci(C2(Grm + Grr)! (31) 


where the coefficient C4 is 0.518 and C2 is 0.61. 

The problem in our work is similar to that described in lit- 
erature [20], therefore, it can be derived that the Sherwood 
number is function of the sum of Grashof number and diffusional 
Grashof number based on above formula and our numerical 
results. 


Sh = 0.65(Grm + Grr)!/4 (32) 


This equation figures out the essential factors that influent the 
oxygen transport, and helpful to understand the interaction of 
heat/mass transfer and dynamics inside the cell. Using the equa- 
tion, the oxygen transport coefficient can be easily calculated, 
and then limiting current density can be obtained very simply. 
In addition, Sherwood number calculated from above equation 
can also be used to evaluate if the configuration of the cell is 
reasonable and to optimize the oxygen transportation in the cell. 

Fig. 11 shows concentration losses caused by oxygen insuf- 
ficiency. The reduction of cell performance due to oxygen 
transport limitation is shown to be proportional to oxygen trans- 
fer coefficient. 


3.4. Effect of water transport 


Generally the water management is vital problem for any fuel 
cells. However, the water transport through the membrane for 
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Fig. 11. Concentration over-potential distribution in catalyst surface. 
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Fig. 12. Distribution of net water molecules per proton flux through the mem- 
brane. 
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Fig. 13. Distribution of ohmic over-potential. 


an air-breathing PEMFC has some particular characteristics in 
comparison with general PEMFC with external humidification. 
The anode side is easily dried without humidification, and the 
sole water source in anode side is water diffusion from cath- 
ode side due to water concentration gradient. Fig. 12 is net 
water molecules per proton flux through the membrane with the 
defined water transport direction from anode to cathode. Without 
external humidity, operating in dry anode gas, a certain amount 
of product water can transports from cathode to anode by back 
diffusion. While protons migrate from the anode to cathode and 
carry with them water molecules when current is drawn from 
the cell. Hence, there is always a little osmotic drag; the amount 
of water collected at the anode is the balance between the total 
amount of water back transported and of water transported again 
to the cathode due to osmotic drag. It is observed that the net 
water flux is negative because of the presence of back diffusion. 
And osmotic drag takes place wherever the current density is 
high because there is more proton transport, and in the other 
regions, back diffusion takes place. It can be seen the osmotic 
drag takes place in the region near the channel opening, and the 
diffusion occurs in the center and anode inlet region. Fig. 13 is 
voltage losses due to membrane proton transport resistance. The 
center region is better hydrated due to stronger water back diffu- 
sion, and the ohmic losses are consequently low. The regions in 
the open area of cathode channel have worse water condition due 
to osmotic drag effect, the region near the anode inlet the ohmic 
losses are high because pure hydrogen is used as anode reactant. 


4. Conclusions 


A complex three-dimensional model including detailed 
heat/mass transfer parameters and phenomena has been devel- 
oped for an air-breathing PEMFC with real geometry and 
structure. The model accounts for all major phenomena except 


phase change. The numerical results for the cell performance 
are physically consistent and be in good agreement with the 
experimental data. The effects of the heat and mass transfer 
under the free convection, and the water transport on the air- 
breathing PEMFC performance have been investigated through 
the examinations of distribution of the velocity, temperature, 
concentration of gas component, dimensionless coefficient for 
heat and mass transfer analysis, current density and the con- 
centration over-potential. The profound interactions between 
velocity, temperature, concentration of gas component and cur- 
rent density are clearly shown from the ‘calculated results and 
their analysis’. The capability of the comprehensive model 
for analyzing and understanding many interacting in an air- 
breathing PEMFC that cannot be studied experimentally is 
demonstrated. As the results this model provides more pro- 
found insights into the interacting phenomena between heat, 
mass transfer and electrochemical reactions in the cell. The 
non-dimensional oxygen transport coefficient (Sh) has been 
calculated by full N-S function solution. And then a simple 
equation of Sherwood number has been derived from the sum of 
the Grashof number and the diffusional Grashof number. This 
method is very useful to simplify the calculation and to optimize 
the design of an air-breathing PEMFC. The model provides a 
computer-aided tool for the design and the optimization of a 
future air-breathing PEMFC with better performance. 
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